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ABSTRACT 

Let  S  denote  the  space  of  bivariate  piecewise  polynomial  functions  of 

degree  <  k  and  smoothness  p  on  the  regular  mesh  generated  by  the  three 

directions  (1/0)/  (0,1),  (1,1).  We  construct  a  basis  for  S  in  terms  of  box 

splines  and  truncated  powers.  This  allows  us  to  determine  the  polynomials 

which  are  locally  contained  in  S  and  to  give  upper  and  lower  bounds  for  the 

degree  of  approximation.  For  p  =»  L( 2k-2 )/3j  ,  k  /£  2(3),  the  case  of  minimal 

degree  k  for  given  smoothness  p,  we  identify  the  elements  of  minimal  support 

in  S  and  give  a  basis  for  S,  »  {f  e  S:  supp  f  C  ft},  with  ft  a  convex 

ioc  ~ 

2 

subset  of  R  . 
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SIGNIFICANCE  AND  EXPLANATION 


Local  support  bases  for  piecewise  polynomial  spaces  are  important  for 
applications  such  as  finite  element  methods,  data  fitting  etc.  In  [BH^ ]  a 
general  construction  principle  for  such  "B-splines"  was  used  to  obtain  the  so 
called  box- splines.  They  have  a  particularly  regular  discontinuity  pattern 
and  coincide  in  special  cases  with  standard  finite  elements. 

— This  report  investigates  the  use  of  translates  of  certain  bivariate  box- 
splines  in  the  construction  of  a  unified  theory  for  piecewise  polynomial 
functions  on  regular  meshes. 

A  simple  mesh  is  considered,  derived  from  a  square  mesh  by  drawing  in  the 
same  diagonal  into  every  square.  The  space  S  of  piecewise  polynomial 
functions  of  a  given  degree  and  smoothness,  and  with  discontinuities  (in  some 
derivative)  only  across  lines  of  that  mesh  is  considered.  He  show  that  the 
box  splines  and  their  translates  provide  a  basis  for  the  ,^local"  part  of  S 
and  use  the  techniques  of  IBH^J  to  analyse  the  approximation  properties  of  S. 

The  report  stresses  the  importance  of  local  support  bases  which  are 
desireable  for  applications  such  as  finite  element  methods,  smoothing  of  data 
and  approximation  in  general.  Our  results  should  be  useful  for  the  further 
investigation  of  smooth  piecewise  polynomials,  in  particular  on  regular  meshes 
(c.f.  [CW]  ,  [Si],  [Si]  for  related  work).  " 


Bivariate  box  •rlines  and  smooth  pp  functions  on  a  three-direction  mesh 

C.  de  Boor  and  K.  HOXlig 

0.  Introduction .  This  paper  records  further  results  of  our  continuing  investigation  of 
certain  multivariate  B-splines.  It  follows  [BH1 ]  in  which  we  discussed  general  properties 
of  box  splines  and  the  spaces  spanned  by  translates  of  a  box  spline. 

In  the  present  paper,  we  explore  the  question  to  what  an  extent  box  splines  nay  be 
useful  in  the  study  of  spaces  of  smooth  pp  (:*  piecewise  polynomial)  functions  in  which 
they  lie.  We  restrict  attention  to  the  simplest  interesting  situation,  that  of  the  space 


of  bivariate  pp  functions  in  CP  ,  of  degree  <  k  ,  on  the  mesh  A  obtained  from  a  uniform 
square  mesh  by  drawing  in  the  same  diagonal  in  each  square.  Even  in  this  simple  setting,  we 
find  much  challenge}  in  fact,  we  must  leave  some  obvious  questions  unanswered  for  the 
present. 

The  specific  questions  we  tried  to  answer  are:  (i)  Are  these  B-splines  "basic",  i.e., 
to  what  an  extent  do  box  splines  provide  a  basis  for  S  ?  The  answer  is  that  they  provide  a 
spanning  set  for  the  "local  part",  but  have  to  be  augmented  by  certain  truncated  powers  to 
give  a  spanning  set  for  all  of  S  .  In  certain  special  circumstances,  they  even  provide  a 
basis  for  all  finitely  supported  elements  of  S  .  But  this  happens  rarely,  because  the 
answer  is  "Usually  not"  to  the  question:  (ii)  Are  these  B-splines  "minimal",  i.e.,  does 
S  contain  no  element  with  support  strictly  inside  that  of  a  box  spline?  The  box  splines  do 
provide  material  help  in  answering  the  question:  (iii)  What  is  the  approximation  order 
from  S  ? 

In  outline,  the  paper  is  as  follows:  In  Section  1,  we  introduce  the  relevant  notation 
in  the  process  of  specializing  the  general  results  of  [BH^ ]  concerning  box  splines  to  the 
specific  context  of  the  bivariate  3-direction  mesh  A  .  We  study  the  space  spanned  by 
certain  translates  of  one  such  box  spline,  prove  these  translates  to  be  linearly 


independent  even  locally,  and  characterize  all  polynomials  in  their  span.  In  Section  2,  we 
show  that  S  is  spanned  by  certain  box  splines  and  their  translates  together  with  certain 
truncated  powers.  These  latter  functions  are  zero  on  a  half space  and  agree  with  a  suitable 
polynomial  on  the  complementary  halfspace.  This  permits  us  to  show,  in  Section  3,  that,  in 
effect,  the  approximation  order  from  s  is  entirely  determined  by  how  well  one  can 
approximate  from  Sloo  :•  span  of  box  splines  contained  in  S  .  This,  in  turn,  can  be 
related  to  the  question  of  which  polynomials  are  contained  in  Sloc  .  We  answer  this 
question  in  full  and  thereby  obtain  upper  and  lower  bounds  on  the  approximation  order 
from  S  which  coincide  in  some  cases  and  are,  in  any  event,  very  close  when  p  is  as 
large  as  possible,  i.e., 

P  -  POO  i-  |j2k-2)/3j  . 

We  also  give  a  conjecture  concerning  the  approximation  order  for  p  <  p(k)  . 

In  section  4,  we  look  for  elements  of  minimal  support  in  S  .  These  are  provided  by 

the  box  splines  in  case  k  =  1(3)  and  p  “  p(k)  .  For  k  =  0(3)  and  p  -  p(k)  ,  there  are 

in  S  elements  of  smaller  support  than  that  of  the  box  splines.  These  were  first  discussed 

by  Frederickson  [Fr]  .  In  either  case,  we  show  that  these  minimal  support  elements  provide 
a  basis  for  all  finitely  supported  functions  in  S  .  We  also  discuss  the  case  k  s  2(3) 
and  p  «  p(k)  in  which  the  degree  is  not  minimal  for  the  given  p  to  illustrate  that  the 
search  for  minimal  support  elements  can  be  quite  frustrating  when  k  is  not  minimal.  Only 
for  sufficiently  large  k  (with  respect  to  p  )  does  the  minimal  support  question  become 
simple  again. 
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1.  Box  splines  ob  •  three-direction  Bosh.  In  [BH,],  the  box  wglia*  M_  is  defined  so 

the  distribution  on  bP  given  by  the  rule 

M _*♦  ( — >  /  ♦(  £  X(i)  c  )  dX  (1) 

a  to.  1]  i-1 

for  some  sequence  =  *=“  <ei)”  ln  *"  •  In  thls  *ection'  "•  specialise  the  general  results 
of  [BH j]  concerning  M_  and  the  span 

S-  «*  *P«n(Ks<— »>)*«,■ 

of  its  integer  translates  to  the  simple  situation 


m  -  2  , 


with  the  three  directions  given  by 


(djir,  d2«s,  dj»t) 


d,  ev  d2  «-  e2,  dj  s-  e,  +  e2  . 

By  this  we  mean  that  ran  =  C  (d,,  d^  d3>  and  that  r,  a,  t  are  the  relevant  direction 

■oltiplicities  which  characterise  a  ,  !■*•« 

r  l{i:  Cj-  d^l  .  a  I U*  Cj-  d2>l  ,  *nd  t  I  U»  ^  d3}|  . 

This  special  choice  of  2  allows  us  to  delve  nore  deeply  into  the  details  in  a  setting  of 
possibly  practical  importance. 

In  later  sections,  we  will  write 

M  _  ,  instead  of  *_  • 

T|f|l  S 

For  the  remainder  of  this  section,  though,  we  write 


and  write 


for  any  particular 


He  now  study 


M  instead  of  M_ 


My  j-  H(*-v) 


v  €  V  * 


S  S_  -  span  (My)v 


S  is  a  subspace  of 


it  :®  pp  functions  of  degree  ^  i  on  the  partition  A  , 

k  ,d 


k  n-2 


3 


and  A  the  partition  of  R*  into  triangles  obtained  from  the  three  families  of  meshlines 

V  +  ,  V  €  V  ,  z  €  R  . 


He  have  foregone  the  opportunity  to  make  the  syimMtries  in  A  more  apparent  by  having  the 

o 


three  families  of  meshlines  intersect  each  other  at  an  angle  of  120  (as  is  done,  e.g.,  in 
(Fr] ).  This  would  needlessly  complicate  the  notation.  It  is  sufficient  to  note  that  any 


permutation  of  the  meshline  families  can  be  accomplished  by  some  linear  map  on  Rz  ,  and 


the  corresponding  change  of  variables  leaves  r.  .  f"\  invariant. 

Kp  A 


The  smoothness  of  M  depends  on  the  direction  multiplicities#  We  have 


m  e  i^cc14-11 


with 


d  -  (n  -  max  {r,s,t}  )  -  1 

the  number  defined  in  [BH^;  (2.6)]  as  evaluated  for  our  special  case.  Since  n  »  r  +  s  +  t 
,  it  follows  that,  for  fixed  degree  k  »  n-2  ,  we  get  marl  me  1,  smoothness  by  choosing 

max  {r,s,t}  «  |"(k+2)/3l  . 

Then,  for  k  ■  3u  +  i  ,  the  corresponding  maximal  d  is 

d(k)  «-  2u  +  i  -  [(2k+1)/3j  ,  (3] 

i— 1,0,1  .  For  k  “  3(1+1  ,  there  is  Just  one  choice, 

r  -  s  ”  t  *  u+1  , 

while,  for  k  «  3|i  or  k  »  3u-1  ,  there  are  three  choices  for  (r,s,t)  . 

Recall  from  [BH^ ;  Cor. 2  of  Thm.5]  that  d  also  governs  which  polynomial  spaces  are 
contained  in  S  •  Precisely, 

»  c  S  iff  m  <  d  .  (4! 

m  ** 

Of  course,  as  we  will  see  shortly,  some  polynomials  of  degree  higher  than  m  may  also  be 
in  S  • 


It  follows  from  (1)  that 

supp  M  -  {  X(i)5i  :  X  €  (0,1]n)  -  {E*  X(i)di  :  X  €  (0,r] x[0,s] *[0,t] }  . 

Thus,  supp  M  is  a  hexagon  composed  of 

N  :»  rs  +  rt  +  st 


triangles  of  A  which  are  translates  of  the  triangle  spanned  by  d1  and  d.  and  a  like 


number  of  triangles  which  are  translates  of  the  triangle  spanned  by  d2  and  d3  .  This 
implies  that  exactly  N  M^s  have  any  particular  triangle  of  A  in  their  support. 

Since 

detfdjy  dj)  »  1  for  i  +  j  , 

we  conclude  from  [BH^Prop.4]  that  (Hy)v  need  not  be  linearly  dependent.  We  now  prove 
much  more. 

Proposition  1.  tHy^yev  is  locally  linearly  independent,  i.e. t 

{Mv  »  supp  &V(\K  ?  0  }  (5) 

is  linearly  independent  over  any  open  set  A  contained  in  some  triangle  of  A  . 

Proof.  Since  (S)  contains  exactly  N  elements,  it  is  necessary  and  sufficient  to 
prove  that  S  contains  N  functions  which  are  linearly  independent  over  A  .  This  latter 
condition  is  shown  to  hold  once  we  show  that  dim  >  N  .  It  then  follows, 

incidentally,  that  dim  (S  As)  ■  N  .  Here, 

SA*  *>  *3 

is  the  linear  space  of  all  polynomials  contained  in  S  . 

The  proof  consists  in  identifying  various  elements  of  *„  .  For  the  specific  S  ,  we 
have  from  [BH ^ (Theorem  5]  that 

w_  -  ker  Aker  t\  ker  •  (6) 

Correspondingly,  we  would  like  to  specify  linearly  independent  elements  of  in  the  form 


f  P<r, 

0<S, 

T-t 

I1Pl2°I3T<1> 

for  (  p<r. 

»*S, 

T<t 

J  ' 

V  P-r, 

P<8, 

T<t  , 

1 

with  ly  ,  12  f  Ij  right  inverses  of  D|  ,  Dj  ,  D^+Dj  respectively,  and  1  denoting  the 
function  x  \ — ►  1  .  But,  since  each  of  these  integral  operators  fails  to  commute  with  at 
least  one  of  these  differential  operators,  it  is  tricky  to  make  the  construction  precise  in 
this  form. 


Instead,  we  single  out  the  two  classes 


A.  >“  (♦«.*  0(1)  <  *■>  and  A,  I-  {*  t  a(  2)  <  ■> 

1  a  t  » 


of  monomials 


♦a  *“  O**/®1  • 

Than  A,  C  ker  D^r  ,  Ajfku  D2S  ,  therefore  a,Ha2  provides  a  linearly  independent  set 
of  re  elements  in  t_  .  In  addition,  we  pick  a  set  Bj  £  span(AjXA2)  of  rt  elements 
and  a  set  B2  C  spanCAjVj^)  of  st  elements  in  «_  and  are  then  certain  of  the  linear 


independence  of  the  total  collection 


(A1f\A2)  WB,  UB2 


os  soon  as  we  prove  that  both  Bj  and  B2  are  themselves  linearly  independent. 

To  construct  B1  ,  we  consider  the  right  inverse  J  of  D,+D2  for  which 

(jf){*,0)  -  0  . 

To  find  J*  ,  we  write  J*  -  Z  ca+0  and  consider  the  resulting  linear  system 

<*  01  P  P 

*  We^-..,’  -  *a 
E  c8#e(.,0)  -  0 

This  gives  c8  *  0  for  0(2)  »  0  and  therefore  c^  “  0  for  |0|  #  |a|  +  1  and  for 
|0|  »  lol  1  with  0(1)  >  a(1)  ,  hence  also  c__  -  1  .  In  conclusion, 

2  ' 


*0  s"  ^0  +  ,p*n  ^  “  181  *  id)  <  8(D)  • 


Therefore,  more  generally. 


J‘V  S  • 


Now  set 


By  (7), 


8,  s-  »  i«1,...,t»  j-0,...,r-1) 


jl*(j,s-1)  €  #<j,.-1+i)  ^  8p*n  ^1NA2^  f°r  1>0'  J<r  ' 

hence  Bj  C  span  (a^\  a2)  ,  as  desired.  This  also  implies  that  B.,  C  ker  D1r  ,  hence 

B^  Cher  Aker  D1r(D1+D2)t  . 


J..-,,)  -  S*«»2’ 


6 


therefore  (D1+D2)t[B1l  C  ker  D2°  ,  hence  also  C  her  D28(D1+D2)t  •  We  conclude  that 


B.  O  e 

1  —  a 


Finally,  we  need  to  show  the  linear  independence  of  Bj  .  For  this#  consider  the 
matrix  C  of  polynomial  coefficients  for  the  elements  of  Bj  #  i.e.# 

jSj  a-1  »x  £  C(i,j»8)  ig  ,  1-1,... tj  j-0,...,r-1  . 

Choose  the  (reverse  lexicographic)  ordering 


'i+j  < 

L,j)  <  (h,k)  s-  1  or 
xi+i  - 


Vi+j  -  h+k  and  i  <  h 

Then  C  is  unit  lover  triangular  in  the  sense  that 

«  +j/S-1+1  ♦  «P»“UB  x  8  <  <j,s-1+i>>  ,  all  j,i, 
hence  of  full  rank.  Thus  Bj  is  linearly  independent. 

The  construction  of  B2  proceeds  in  exactly  the  sane  vay,  with  the  roles  of  the  two 
independent  variables  interchanged.  1 | | 


i,  take  (r,s,t)  -  (2,2,1)  .  Then  n-r  +  s  +  t-  5  ,  hence  k  -  3  .  Also, 


d  -  d(3)  -  2  ,  hence  M  is  a  piecewise  cubic  C  function.  How 


A^Aj  -  x  a(1),  a( 2)  -  0.1  } 


forms  a  basis  for  *  ,  x-  bilinear  polynomials.  Further, 


Bj  -  x  i-li  j-0,1  }  , 


J+0,1  “  +0,2  J+1,1  "  +1,2  +  C+0,3 


We  determine  c  from  the  condition  that  (D1+Dj)J^1  (  ^  (  .  This  gives 


hence  c  -  -1  .  Thus  B 


+0,2  +  +1,1  +  C+0,2  ”  +1,1  ' 

1  *  {+0, 2  '  +1,2  ’  +0,3}  *  ■*  try' 


B2  "  W2  „  '  +2,1  “  +3,0^  *  Th*r«for*»  Altogether, 


*S  "  *2  +  ,pan^  +  1,2  “  +0,3  '  +2,1  "  +3,0* 


in  case  S  -  (e^,  et,  e2#  e2#  ejtej)  . 


Corollary  1.  For  any  triangle  t  of  A  ,  a  baala  for  »_  la  provided  by  the  N 
nontrivial  polvnomiala  which  agree  with  on  t  . 


Corollary  2.  _I£  pen  agrees  with  M  on  some  triangle  of  A  ,  then  p  e  . 

H— ark.  While  Corollary  1  is  quite  special.  Corollary  2  is  valid  for  an  arbitrary  box 
spline  M  in  any  number  of  dimensions.  This  is  a  consequence  of  [BH^ ;  Theorem  S]  and  is 
due  to  the  fact  that  any  polynomial  p  which  agrees  with  M  on  some  open  set  is 
necessarily  mapped  to  0  by  any  differential  operator  Dz  for  which  DZM  is  supported 
only  on  certain  hyperplanes. 


2.  Spanning  nti  and  local  bases.  In  this  section,  we  give  a  truncated  power  basis 


for  (Q)  ,  with  Q  any  rectangle  bounded  by  A-mesh  lines.  This  would  allow  us  to 

verify  the  dimension  formulas  of  [CW1 ]  for  this  space.  He  also  give  a  spanning  set  for 
itself  which,  though  finitely  linearly  independent,  permits  nontrivial  infinite 

k,a 

linear  combinations  which  add  to  zero.  Its  main  feature  is  that  it  consists  of  finitely 
supported  functions,  viz.  box  splines,  on  the  one  hand  and  of  functions  supported  on  half 
spaces  and  agreeing  with  some  polynomial  there  on  the  other.  These  latter  functions  do  not 
contribute  to  the  approximation  order  obtainable  from  the  scale  ,  as  we  will  show 

in  the  next  section.  This  means  that  the  approximation  order  is  no  better  than  that 
obtainable  from  the  span  of  the  relevant  box  splines,  and  this  fact  allows  us  to  give  upper 
and  lower  bounds  on  the  approximation  order  which  differ  by  at  most  two  in  case  p  is  as 
large  as  possible. 

Recall  from  [BH^]  th_» c  Dahmen's  truncated  powers  [D]  can  be  thought  of  as  shadows  of 

the  standard  cone  :  With  S  (5^)"  in  R?1  ,  the  corresponding  truncated  power  or 

cone  spline  C_  is,  by  definition,  the  distribution  on  given  by  the  rule 

2 


c3  i  ♦  !— >  /  n  ♦(  E"  MiKj  ax 


Since  M-:$  | — >  /  $(l™A(i)£  )  dX  ,  it  follows  that 

=  [0,1]"  1  1 

C_ 


l  «_(•  -  v) 
vez^  a 


(i> 


Recall  from  [BHj ]  that,  for  Z  in  a  , 


DC-  «  C_.  „  . 
Z  s  s\Z 


(2) 


Now  specialize  to  the  setup  of  Section  1,  i.e.,  to  the  specific  sequence 


(d1 sr,  d2ss,  d3«t) 


consisting  only  of  the  vectors  d1  ■  e1 ,  d2  =  e2  ,  and  dj  =  ei+e2  in  *  »  and  therefore 
characterized  by  the  corresponding  direction  multiplicities  (r,s,t)  .  It  follows  from  (2) 
that 


C-  has  all  derivatives  of  order 


|s+t-2| 
<  <  r+t-2  l 


continuous  across  span! 


k\>  • 


(3) 


Correspondingly,  the  univariate  function  N_  given  by  the  rule 

N„(z)  t«  C_(z,1-z)  ,  all  z  «  K  , 
s  s 

ia  a  univariate  B-spline,  i.e., 

C_(z,1-z)  “  c._M(z|0jr,  V2  :t,  lie)  (4) 

for  eoee  positive  c_  . 

Here  is  an  outline  of  what  ia  to  follow.  He  show  that,  near  a  lower  left  corner  of  its 

support,  any  f  €  wj*  can  be  written  as  a  linear  coabination  of  certain  truncated  powers. 
k»o 

For  this,  we  split  f  into  its  homogeneous  components.  Being  homogeneous,  each  such 
component  is  determined  by  its  restriction  to  a  line  which  "cuts  across"  the  corner.  Such  a 
restriction  is  a  univariate  spline,  hence  uniquely  representable  as  a  linear  combination  of 
certain  univariate  B-splines,  i.e.,  of  restrictions  of  certain  truncated  powers. 

Next,  on  subtracting  from  f  this  linear  coabination  of  truncated  powers,  we  obtain  a 
0 

new  element  of  z  .  whose  support  is  inside  that  of  f  and  offers  lower  left  corners,  to 
K,  A 

the  right  and/or  above,  for  further  "peeling  off". 

He  begin  with  a  study  of  the  simple  pp  apace  which  models  the  behavior  of  f  €  A 
near  a  lower  left  corner  of  its  support.  He  denote  this  space  by 

S(<k,v)  I* 

2 

and  mean  by  this  the  space  of  all  pp  functions  of  degree  <  k  with  support  in  R+  and 
possible  singularities  only  across  the  three  rays 

*tdi  ,  i-1,2,3  . 

In  addition,  we  think  of  v  here  as  a  3-vector,  with  v(i)  indicating  that  all 
derivatives  of  order  <  v(i)  are  required  to  be  continuous  across  B+d^  ,  i»1, 2,3. 

Let 

Ht  i-  {  f  s  f ( zx)  -  zlt(x)  ,  all  x  €  R2,  z  €  R+  > 
denote  the  collection  of  all  functions  on  R2  (positively)  homogeneous  of  degree  t  .  As  is 
well  known, 

z.  -  9  (H.nO  , 
t<k 

so  it  makes  sense  to  talk  about  the  homogeneous  component  of  degree  t  of  a  polynomial.  He 


now  ukt  the  same  claim  for  S(<k,v>  .  With 


we  claim 


S(i,v)  Ht/^S(«k,v)  , 


S(<k, v) 


-  •  S(t,v)  . 

Kk 


Proof.  Only  tha  inclusion  requires  proof.  To  prove  this  inclusion,  it  is 
sufficient  to  show  that  the  pp  function  made  up  of  the  ith  degree  homogeneous  components  of 
an  element  of  S(<k,v)  is  again  in  s(<k,v)  .  This  follows  from  the  following 


Claim.  If  a  polynomial  p  vanishes  to  order  p  along  the  ray  Rj.d  ,  l.e., 

D  p  •  0  on  R+d  for  [  a|  <  p  , 

then  each  of  its  homogeneous  components  also  vanishes  to  order  p  on  K*  • 


Proof.  Assume  without  loss  that  d  *  e2  .  Since  p  ■  0  on  R^ej  for  l  <  p  ,  we 
must  have 


,.{p+1,0> 

p  "  O  q 


for  some  polynomial  q  .  Writing 


p  ■*!  Ip.  with  p,  e  H.  ,  all  i  , 

(<|r  *  *  * 

we  conclude  that  each  has  the  factor  {)'’”'  ,  therefore  vanishes  on  Re2  together 

with  every  derivative  of  order  <  p+1  .  (In  particular,  p^  -  0  for  t  <  p+1  ).  ||| 


We  took  the  trouble  to  express  S(<k,v)  in  terms  of  its  homogeneous  components 
S(t,v)  since,  on  S(i,v)  ,  the  linear  map  R  given  by  the  rule 

(Rf)(x)  f(z,1-z)  ,  all  z  e  R  , 

is  1-1.  This  follows  from  the  fact  that,  for  any  f  e  , 

f(Xz,X(1-z) )  -  X*(Rf)(s)  ,  all  X  e  R+  ,  z  8  R  , 

hence  such  f  is  determined  on  the  entire  halfspace  x(1)  +  x(2)  >  0  once  Rf  is  known. 
We  claim  that  R  carries  S(t,v)  onto  the  univariate  spline  space 

8u(<t,v) 

1 1 


which  consists  of  all  pp  functions  g  of  degree  <  £  on  R  with  breakpoints  0>  1/2,  1, 


with  support  in  [0,1]  ,  and  with 


v(i) 

gee  near  5 


‘■-{’8  “  *-0 


Indeed,  R  carries  all  of  S(<£,v)  into  .  In  addition,  we  recall  from  (4)  that. 


3  «  (d1  jr,  djts,  Sjit)  , 

R  carries  the  cone  spline  C_  to  a  positive  multiple  of  the  univariate  B-spline 


M(*|0tr,  1/2st,  Its) 


This  implies  that,  with 


(C±)  «-  (d^E-vd),  d3»l-v<3),  djt  l-v(2) )  , 


the  cone  splines  C 


are  in  S(£,v)  ,  and  R  carries  these  to  a  basis  for 


Su<<£,v)  .  Consequently,  these  cone  splines  must  form  a  basis  for  s( £, v)  .  in  particular. 


dim  S<£,v)  -  (E,  <£-v<i>>+  -  £-1  )+  . 


Therefore 


dim  s(<k,v)  -  Et<J{  (E^  (£-v(i>>+  -  £-1  )+  •  (5) 

Raster  k.  This  formula  shows  that  ^  contains  no  finitely  supported  functions  unless 

p  is  suitably  small:  If  f  e  if  A  0  has  finite  support,  then  its  support  must  contain 

K,a 

a  "lower  left  corner",  i.e.,  a  mesh  square  i-  [v(1 ),v(1 )+1 ] »[v(2) ,v(2)+1 ]  along  whose 
left  and  lower  edge  f  vanishes  to  order  p  .  This  implies  that  f  agrees  on  Qv  with 
g(  *-v)  for  same  g  e  s(<k,p,p,p)  .  This  in  turn  implies  that  dim  s(<k,p, p,p)  >  0  ,  and, 
by  (5),  this  is  equivalent  to  having  p  <  (2k-2)/3  .  This  conclusion  was  reached  in  [BD] , 
using  the  same  simple  argument  of  cutting  across  such  a  lower  left  corner  of  the  support, 
as  is  used  here,  we  realized  only  recently  that  this  conclusion  can  already  be  found  in 


We  are  now  ready  to  give  a  cone  spline  basis  for 

*JU(Q)  \,d  |q^cP(«>  ' 

We  use  the  translation  map  given  by  the  rule 

( T  f ) (x)  i»  f ( x  -  v)  . 
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Proposition  2.  Let  Q  “  [0,M+1 1  *(0,N+1  ]  .  Then  ^(Q)  is  the  direct  sun  of  the 


ssaces  V^'V  |Q  w— 


(-1,-1, p)  , 

V  - 

(0,0) 

(-1.P.P)  , 

V  - 

(x,0 ) 

(p*-1,p)  . 

V  - 

(0,y) 

(p»p,p>  , 

V  * 

(x,y) 

,  all  v  e  VQ  s-  V  A([0,M]  x[0,N)  )  ,  and  x,y  >  0 


Consequently 


dim  «£  ^(Q)  -  EUk  i+1  +  (l-p)+  +  2(M+N)  ( t-p)+  +  MN(3(*-p)+  -  i-1  }+  , 
and  the  restriction  to  Q  of  the  cone  splines 

Cr  _  (-v)  ,  all  i,  with  (c  )  (d  :*-v  (1),d  :l-v  (3),d_:t-V  (2)) 

V'^i+i+l  3 


for  t  <  k  ,  v  e  VQ  , 


forms  a  basis  for  wp  . (0)  . 
-  k,A 


Proof.  For  any  choice  of  ,  the  spaces  T^s(<k, v^) ^  are  in  direct  sum.  For  the 

A 

specific  choices  of  given,  they  are  all  In  ^(Q)  .  Thus  it  only  remains  to  show 

that  sp  a<Q)  is  contained  in  their  sum. 

We  proceed  by  induction.  For  this,  we  use  again  the  (reverse  lexicographic)  ordering 

(I v|  <  |w| 
or 

I v|  -  |w|  and  v( 1 )  <  w(1 ) 

which  provides  a  total  ordering  for  Vg  .  We  again  use  Qv  to  denote  the  unit  mesh  square 

whose  lower  left  corner  is  v  .  The  induction  hypothesis  to  be  advanced  is  the  following: 

For  all  v  <  w  ,  there  exists  f„  e  t  S(<k.v  )._  so  that 
-  -  v  v  v  ia  - 

6  f  -  E  f 

W  V<W  V 

vanishes  on  Qv  .  In  order  to  advance  this  hypothesis,  we  now  show  that,  in  its 

v<w 

consequence,  6  agrees  on  Q,  with  some  f  6  T  S(<k,v  )  .  There  are  four  cases: 

(i)  w  «  0  .  Then  5^  »  f  ,  hence  it  agrees  with  some  fw  e  S( <k,-1 ,-1 , p)  . 

(ii)  w  -  (i,0)  for  some  i  >  0  .  Then  5  vanishes  to  order  p  on  the  segment 


1  3 


ix[0,1]  ,  therefore  agrees  on  Q,,  with  some  fw  e  T^S^k.-l  ,p,p)  . 

(Hi)  w  -  (0,j)  for  som  j  >  0  .  Then  6^  vanishes  to  order  p  on  the  se^ent 
[0,1]  «j  ,  therefore  agrees  on  Q*  with  som  fw  e  T^St  <k,  p,-1  ,  p)  . 

(iv)  w  “  (i,j)  for  i,J  >  0  .  Then  6 vanishes  to  order  p  on  the  left  as  well  as 
on  the  lower  boundary  salient  of  ,  therefore  agrees  on  with  som  fw  6 

TwS(<k,p,p,p)  . 

Since  supp  fw  Q  t  r+2  g\  U  fiy  >  this  advances  the  induction  hypothesis,  since  it 

v<w 

implies  that  6  -  f  vanishes  also  on  as  well  as  on  \J  • 

**  »  v<w 

The  diMnsion  formula  (6)  now  follows  from  (5)  .  Ill 


Since  M  and  N  are  arbitrary  positive  integers,  we  obtein  the  following 


Corollary.  For  0  »  R^.2,  .(Q)  is  spanned  by  the  restriction  to  Q  of  the  cone 


splines  listed  in  (7)  (with  V--^2). 


Next,  we  investigate  the  relationship  of  *£^((3)  to  ^  1^  .  These  two  specea  are, 
in  fact,  the  ssm,  but  this  is  not  clear  a  priori.  It  is  obvious  that 

Wk,A|g  €  VaIQ^*®  Wk,A(Q)  • 

For  the  converse  containMnt,  it  is  necessary  and  sufficient  to  show  that  every  f  C 

»?  *(Q>  can  be  extended  to  an  eleMnt  of  if  .  .  By  Proposition  2,  this  is  established 
kf  A  kf  a 

once  we  show  that,  for  each  v  e  Vg  ,  TyS( <k, vv> |q  can  be  extended  to  a  subset  of  , 

and  this  is  obvious  as  long  as  *  (p,p,p)  .  This  leaves  three  cases: 

(i)  vv  "  (— 1 , p ,p )  .  Then,  for  i>1,>..,pt1, 

C  C 

i  V***'CH-t+1 


involves  the  direction  d^  i+2-i  times,  l.e.,  more  than  t-p  times,  hence  fails  to  be  in 
CP  across  Rj.dj  .  Recall  from  (4)  that  the  restriction  RC^  of  given  by 

(Rci)(t)  i-  ct(t,i-t)  ,  all  ten, 

is  a  scaled  univariate  B-spline  involving  just  the  three  knots  0,  1/2,  and  1,  and  the 
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latter  two  no  more  than  t-p  tinea .  We  can  therefore  write  RC^  on  [0,1]  aa  a  linear 
combination  of  the  truncated  powera 

(1/2  -  •)/  ,  <1  -  •)/  ,  r-p*1,...,t  . 

2 

Since  CA  is  homogeneous  of  degree  1  ,  this  implies  that,  on  R^,  itself  ia  a 

linear  combination  of  the  truncated  powers 

Ti,t,rl  ^dix)t"r(dix)+r  *  (8) 

i  “  2,3,  and  r  “  p+1,...,i  ,  with 


d2  *i*  d3  *“  *re2 ' 


and  these  truncated  powers  are  all  in  ^  •  We  conclude  that 

on  *+2  ,  S(<h,wv)  £  S( <k,p,p,p)  ♦  spant^  ^  i  i-2,3»  p<r<t<k> 


(9) 


(ii)  -  (p,-1,p)  .  In  this  case,  we  conclude  that  we  can  write  the  offending  cone 
splines  aa  linear  combinations  of  the  truncated  powers 

*f->  (d^)4-*^)/  and  x|— >  (djxJ^C-djx)^, 

p<r<t<k,  with 

d1  *2  ‘ 


This  implies  that, 

on  *+2  ,  S(<k,vy)  C  S(<k,p,p,p)  ♦  span{T1>t^r«  i-1,3»  p<r«i<k)  +  ,  (10) 

the  last  summand  because  the  function  x  I— >  (d,x)4  r(-d,x)  r  is  in  span{T,  ,  }  ♦  w  . 

(iii)  -  (-1,-1, p)  .  For  this  case, 

S(t,wy)  -  Htn*  on  *+2 

for  t  <  p  .  For  t  >  p  ,  has  either  d^  or  dj  but  never  both  appearing  more  than 

t-p  times.  This  implies  that 

S(<k,vv)  C  »k  +  S(<k,p,p,p)  +  S((k,-1,p,p)  +  S( <k,p,-1 ,p)  . 

Thus,  using  the  other  cases,  we  find  that, 

on  R+2  ,  S(<k,vv)  C  S(<k,p,p,p)  ♦  spanCTj^  ^  »  i-1,2,3>  p<r<t<k)  ♦  .  (11) 


This  establishes  most  of 


(12) 

(13) 

(14) 

(15) 


Proof.  We  only  need  to  prove  (13).  But  this  follow*  from  (12),  tram  Proposition  2,  and 
from  (9)— (11),  using  (1)  to  convert  the  cone  splines  Into  linear  combinations  of 
corresponding  box  splines.  I | | 


It  is  easy  (but  perhaps  not  all  that  useful)  to  obtain  from  Theorem  2  a  spanning  set 

of  the  same  form  for  all  of  if!  .  .  let  ft, ,  ft  ,  ft  be  the  three  domains  into  which  R2 

2 

is  subdivided  by  the  three  rays  ^d^ ,  ^d2,  end  *-d3  ,  with  ft^  “  R+  .  For  given  f  e 

Wk  A  '  w*  may  cho°8e  Theorem  2  an  f^ew^  +  T+  8  which  agrees  with  f  on  (^  .  The 

function  f  -  f,  is  in  w?  .  and  vanishes  on  0.  ,  hence  vanishes  to  order  p  on  30,  .  In 
•  k,A  1  1 

effect,  the  component  ot  f  i  from  ♦  T  insures  that  f  -  fljT  vanishes  to 

order  p  on  3Q  .  This  makes  it  possible  to  write  f  -  f,  -  on  ft.  as  an  element  f ,  _ 

1  »#*1  I/O 

of  the  span  of  the  box  splines  listed  in  (15).  An  analogous  argument  therefore  establishes 
the  existence  of  an  element  f 2(T  in 

spent  x| — >  (di(x-*_di>)*  r(-djx)+r  *  i“1 »  2i  p<r<l<k  ) 

£  ♦  spent  TA  4  (•  -  S_d^)  :  i-1,2j  p<r<i<k  }  £  +  T 

so  that  g  f  -  f 1 >T  -  f2fx  vanishes  to  order  p  on  *-d3  •  Since  f2,T  vanishes  on 
ft,  ,  it  follows  that  g  vanishes  to  order  p  on  3ft,  and  3ft,  ,  while  g  -  f  -  f,  _  on 

1  2  3  1/* 


Qj  ,  hence  9  *  f  1  j  there.  This  ukai  it  possible  to  write  g  on  0^  as  a  linear 
combination  fi<s  of  box  splines  whose  support  is  entirely  in  and  which  are  obtained 

from  the  box  splines  listed  in  (15)  by  the  linear  change  of  variables  which  carries 
to  Q  .  Because  we  chose  the  direction  dj  to  appear  exactly  i-p  times  in  the  box 
splines  of  degree  l  in  (15),  the  box  splines  on  so  obtained  are,  in  general,  not 
translates  of  the  ones  in  (15),  but  could  be  written  as  infinite  linear  combinations  of 
such  translates. 

It  follows  that  f  -  f1(J  +  f2  T  *  f1>s  ♦  f2  g  +  f 3  s  .  In  this  fashion,  we  can 

A 

represent  i'  ^  as  the  span  of  ,  the  truncated  powers  listed  in  (14),  and  certain  box 
splines. 


3.  ifprniaitlaa  ardtr.  in  this  sac  t  ion,  we  give  up  par  and  lower  bounds  for  tha 
approximation  ordsr  of  tha  smooth  pp  spaca 

S  *k,A  ' 

with  A  tha  thraa- dir action  mash  introducad  in  Saction  1.  Tha  approximation  order  of  S 
is,  by  definition,  tha  integer  m  for  which  the  following  holds:  For  all  sufficiently 
smooth  functions  f  , 

dist(f,  Sjj)  ■  0(h“) 

m 

whiles  for  »oaa  C  -function  f  , 

disttf,  S,,)  *  o(h*)  . 

Bara,  tha  scale  (Sh)  of  approximating  spaces  is  generated  from  s  by  simple  scaling, 

®h  *-  vs>  ' 

with 

(ohf)(x)  i*  f(x/h>  ,  all  f,  x,  h  . 

Further, 


dist(f,  0)  i“  infueo  If  -  ul  , 
and  I •  I  is  tha  sup  norm  on  some  closed  domain  OC  I2  , 

If  I  i-  ,m  ““Pago  lt<*H  • 

In  this  definition,  tha  approxiawtion  order  depends  on  8  ,  and  rightly  so.  If,  e.g., 

all  tha  elements  of  S  had  their  support  in  ,  than  would  be  entirely  unable  to 

2 

approximate  to  functions  having  some  support  in  1^  ,  hence  might  well  have  different  - 

approximation  order  depending  on  whether  or  not  8  lies  entirely  in  •  for  the 

specific  spaces  x£  .  or  S_  of  interest  here,  though,  the  approximation  order  is  the 

KpO  2. 

same  for  any  closed  and  bounded  8  with  some  interior,  since,  for  sufficiently  small  h 
,  8^  is  invariant  under  a  suitable  linear  change  of  variables  carrying  one  such  8  into 

another . 

Here  is  a  simple  necessary  condition  for  the  approximation  order  to  be  m  . 


3.  Let  0  be  a  locally  compact  linear  space  of  functions  on  rf1  ,  let  0  be  a 


closed 


let  of  having  0  in  its  interior,  and  assume  that 


diat(f,  ty  -  0(h“) 


This  simple  necessary  condition  la  far  from  aufficiant,  obviously •  For  example, 

taking  a  *  w  ,  wo  obtain  Uj,  ■  0  ,  all  h  ,  and  thia  scale  haa  approximation  order  0  . 

It  ia  not  clear  at  preaent  what  other  conditions  one  should  add  to  get  necessary  and 

sufficient  conditions  for  the  approximation  order  to  be  at  leant  a  .  Yet,  Lemma  3  in 

conjunction  with  Theorem  2  leads  to  a  close-to-exact  estimate  of  the  approximation  order 

of  S  -  w?  .  in  case  p  is  maximal  for  the  given  k  • 
k,A 

A  lower  bound  for  the  approximation  order  of  S  can  already  be  found  in  [BDj  Theorem 
4]  where  it  is  shorn  that,  for  all  sufficiently  smooth  f  , 

dist(f ,  Sjj)  -  0(hP+2)  in  case  p  -  p(k)  [J2k-2)/3j  .  (1) 

We  saw  already  in  Section  1  that  this  is  aa  large  a  p  (•  d-1 )  as  we  can  choose  and  still 
have  box  splines  in  S  •  Correspondingly,  it  is  shown  in  (BD]  that,  for  p  >  p(k)  ,  S  has 
approximation  order  0  . 

Obviously,  (1)  provides  also  the  lower  bound  p(k)+2  for  the  approximation  order  of 
8  in  case  p  <  p ( k )  .  But  one  would  expect  the  approximation  order  of  S  to  increase  aa 
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p  decreases.  This  increase  cannot  be  seen  merely  by  studying  the  approximation  order  of 
S.  Cs  (the  way  (1)  is  obtained).  For,  [BH^Sec.6]  shows  that  the  approximation  order  of 
S_  span(M-( •-X2) )  is  d+1  ,  with 

d  -  k  +  1  -  max(r,s,t) 


in  case 


S  -  (djjr,  djts,  d3it)  and  r+s+t  »  k+2  , 

while,  as  discussed  in  Section  1,  S_  C  ^  and  no  better.  This  swans  that  a  decrease 

in  p  increases  the  number  of  different  box  splines  M_  in  S  ,  but  the  approximation 

order  of  the  additional  subspaces  S.  is  less  than  the  approximation  order  of  those 

a 

already  in  s  when  p  is  maximal,  i.e.,  when  p  -  p(k)  . 

It  is  the  main  goal  of  this  section  to  provide  an  upper  bound  for  the  approximation 
order  of  S  which,  in  the  case  p  -  p(k)  ,  is  close  to  the  lower  bound  (1).  The  proof  idea 
is  simple:  We  show  that  the  polynomials  and  truncated  powers  in  S  do  not  contribute  to 
its  approximation  order.  This  means  that  the  approximation  order  is  already  determined  by 
the  span 

8loc 

of  the  box  splines  in  S  ,  and  the  discussion  in  Section  1  of  the  spaces  S_  spanned  by 

a 

the  translates  of  one  such  box  spline  M_  allows  us  to  determine  the  maximal  m  for  which 
W(b_1  is  contained  in  S^oc  . 

In  view  of  (2.15),  we  set 


with 


S,  I.  S_ 

loc  Kk  i  a 


i,* 


Si,t  *“  ('**'•  ,,,5i+t+l)  and  *“  d3,t_p'  d2,t_(  * 


(2a) 


(2b) 


Proposition  3.1.  I£  dist(f,  3^)  -  0(h*)  for  all  sufficiently  smooth  functions  f  , 

then  w  ,  C  S,„  . 

-  m-1  —  loc 


Proof.  Let  p  be  a  sufficiently  smooth  function.  By  assumption,  there  exist  const 


and  sh  e  Sh  ,  all  h  ,  so  that 


I 


I 


By  Theorem  2,  wa  can  write 


Ip  -  s^l  <  const  h  . 

•h  "  *h  +  "h 


with 


“d  “heV%  +  T) 

and  T  the  span  of  certain  truncated  powers, 

T  ■  span{  T.  .  ( •  -  Xd. )  s  i-1,2,3;  p<r<t<k  >  , 


where 


Ti,i,r<x)  (dix^r^diXlr  * 

Therefore,  for  any  positive  n  ,  the  linear  map 


n  ^  udt  nd2  nd3J 


with 


Ayf  I-  f (  *+y)  -  f 


carries  all  of  o.T  to  zero  (since  k-p  <  k+1 ) ,  as  well  as  all  of  w.  .  He  conclude  that 

n  K 


'V*  W  <  23<lt+1>Con.t  h"  , 


with  O'  »"  (x  e  0  s  diat(x,  38)  >  3(k-p)n  )  . 


For  h  e  n/U  ,  we  have  A  s'  e  o.(S,  )  .  Xn  view  of  Leona  3,  it  therefore  suffices  to 

n  a  IOC 


show  that 


i  ,  C.  ran  A  . 
m-1  -  n 


But  this  is  obvious  since,  for  any  y  f  0  and  any  r  ,  A^  maps  onto  . 


Theorem  3.  For  p  <  p(k)  ^j2k-2)/3j  ,  the  approximation  order  m  of  S  ■  ^ 

satisfies  m  e  [p(k)+2,  m(k)]  ,  with 

m(k)  i-  min {2 (k-p) ,  k+1)  . 


Proof.  The  lower  bound  for  m  is  provided  by  (1).  By  Proposition  3.1,  the  upper  bound 
is  established  once  we  show  that 


m(k) 


«£  S 


loc 
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Obviously,  *  /  8loc  .  further,  by  (2),  each  box  spline  in  sloc  involves  each  of  the 

carries  f  € 

8jqc  to  a  distribution  supported  only  on  the  meshlines  of  A  .  Consequently,  cannot 

contain  the  particular  polynomial  ( ) ^  which  is  carried  by  (d^Dj)*  ^  to  the 

nontrivial  constant  function  x  \~>  (k-p) l2  .HI 

The  next  proposition  shows  that  the  upper  bound  in  Theorem  3  is  sharp  in  the  sense  of 
Proposition  3.1. 

Proposition  3.2.  For  sloc  as  given  1st  S2L. 

max{m  :  «  ,  C  S,  }  -  m(k)  min{2(k-p),  k+1 }  .  (3) 

—  loc 

Proof.  The  assertion  is  obvious  for  k  -  0,1,2.  Since  we  already  know  that 

*m(k)  *  Sloc'  W®  only  httV*  40  “how  that  *B(k)-1 ^  Sloc  * 

Consider  the  box  splines  listed  in  (2).  Specifically,  choose  i  so  that  3^  5^ 

involves  direction  d.,  as  little  as  possible  yet  at  least  once.  Let  (r,s,t)  be  its 

direction  multiplicities.  Then  r  >  1  .  Further,  s+t  -  2(k-p)  as  long  as  k+2  >  2<k-p)  . 

In  the  contrary  case,  r  ■  1  and  s+t  -  k+1  .  Hence,  in  either  case, 

r  -  k-p  ♦  1  -  i  >  1  and  s+t  “  m(k)  .  (4) 

We  conclude  with  (1.6)  that 

()°  e  S-  for  |a|  <  m(k)  and  a(1)  <  r  . 

Indeed,  each  such  polynomial  is  in  ker  ,  and,  as  a  polynomial  of  degree  <  m(k)  “  s+t  , 
trivially  in  ker  D2*<D1+D2)t  * 

Assume  that  we  already  know  that 

()a  e  S1qc  for  )a|  <  m(k)  and  a(1)  <  r+j  (5) 

(as  we  do  now  for  j  ”  0).  if  k-p  >  r+j  ,  then,  with  (4),  k-p  >  k-p  +  1  -  i  +  j  ,  or, 
i-j  >  1  ,  hence  we  may  consider  •  bet  (r+j+1,  s',  t' )  be  its  direction 

multiplicities.  Then 

s'+t*  -  m(k)  -  j  -  1  . 
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By  th«  proof  of  Proposition  1 ,  s_  contains  an  element  of 

s 

»■  0/61  +  span{  ()'  »  |y|  -  |  B| ,  y(1)  <  6(1)>  (6) 

provided  <(1)  <  r+j+1  and  6(2)  <  s*+t'  -  m(k)  -  j  -  1  .  This  implies  that  S_ 

“i-j-1 

contains  an  element  of  (6)  for  B(1)  *  r+j  and  | 61  <  m(k)  ,  and,  on  combining  this  with 
(5),  we  conclude  that  (5)  holds  also  with  j  replaced  by  j+1  . 

This  allows  us  to  conclude  that  (5)  holds  for  r+j  -  k-p  .  But  then,  by  symmetry,  also 
()°  e  S1qc  for  |a|  <  m(k)  and  a(1 )  >  k-p  , 
and  this  finishes  the  proof.  | | | 

For  p  -  p(k)  ,  the  bounds  in  Theorem  3  are  particularly  tight.  If  i  ■  -1,  0,  1 
and  k  ~  3p  +  i  ,  then  p(k)  —  2p-l  +  i  .  Therefore,  for  k  •  3u  +  1  ,  the  approximation 
order  of  »£***  equals  p(k)+2  -  2y+2  .  For  k  -  3u  ,  it  lies  between  p(k)+2  and 
p(k)+3  .  The  particular  case  k  »  3  is  discussed  in  detail  in  (BH2)  where  it  is  shown 
that  the  approximation  order  of  &  is. only  3  («  p(3)+2)  rather  than  4  (»  m(3)). 

This  is  surprising  and  disappointing,  since  it  shows  that  the  simple  mechanism  on  which 
Proposition  3.1  is  based  is  not  sufficient  to  predict  the  approximation  order.  One  might  be 
tempted  to  conclude  from  this  example  and  from  the  case  k  -  3p+1  that  the  highest 
approximation  order  obtainable  from  an  S_  in  S  determines  the  approximation  order  of 
S  itself,  at  least  when  p  -  p(k)  .  The  simple  example  k  -  2  contradicts  this.  In  this 

case,  p(k)  «  0  ,  i.e.,  p(k)+2  »  2  ,  yet  local  polynomial  interpolation  is  well  known  to 

provide  approximation  order  3  from  continuous  piecewise  parabolic  functions  on  any 
triangulation  A  . 

For  p  <  p(k)  ,  the  lower  bound  stays  constant  while  the  upper  bound  increases  until 
it  reaches  k+1  ,  exactly  at  the  point  where  p  is  small  enough  so  that,  already  for  a 
two-direction  mesh  Z  ,  j.  contains  finitely  supported  functions  (see  [BD]).  We 
conjecture  that  the  approximation  order  of  ,  never  differs  from  its  upper  bound 
m(k)  by  more  than  1  .  Proof  of  this  conjecture  would  require  construction  of  a  local 

approximation  scheme  which  makes  use  of  much  of  Sloo  rather  than  just  one  S_  . 
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4.  minimality  of  rapport,  in  this  Motion,  we  show  that  box  splines  in 


Bay  or  nay  not  have  minimal  support,  even  in  the  very  restricted  setting  of  maximal 
smoothness ,  i.e.,  when 

P  -  P(lc)  , 

as  we  assume  throughout  this  section.  Precisely,  we  show  that,  for  k  £  1(3)  ,  the  sole  box 
spline  in  S  has  minimal  support,  while,  for  k  £  0(3)  ,  the  box  splines  in  S  do  not.  In 
the  latter  case,  we  show  that  a  certain  element  first  constructed  in  (Fr]  has  minimal 
support  as  does  its  'flip*.  In  either  case,  the  minimality  allows  us  to  conclude  that 
translates  of  the  minimal  support  elements  span  the  subspace  of  all  finitely  supported 
elements  of  s  .  For  the  special  cases  k  -  3  or  4  ,  this  has  also  been  proved  in  [CW21 . 

The  final  case,  k  £  2(3)  ,  gives  a  hint  of  the  complications  awaiting  those  wishing  to 
study  the  minimal  support  question  for  arbitrary  p  .  We  merely  discuss  the  specific 
choices  k  «  2,  5,  8,  and  state  some  conjectures  concerning  the  general  cara. 

We  make  use  of  the  notation 

for  the  box  spline  M_  with  3  ■  (d^  ir,  d^s,  d^it)  . 

We  say  that  f  has  minimal  support  in  8  if  ft!  and  tSm-  «nly  g  «  X  having 
support  strictly  inside  supp  f  is  g  “  0  .We  My  that  f  has  unique  minimal  support 
in  8  if  f  e  S  and  any  g  e  S  having  support  in  supp  f  is  a  multiple  of  f  .  Clearly, 
any  f  having  unique  minimal  support  in  S  has  minimal  support  in  S  . 

Theorem  4.1.  Let  k  -  3p-2  and  p  ■  p(k)  (-  2p-2)  ,  and  set  S1  ^  .  then  the 

box  spline  M  s-  ^  y  has  unique  minimal  support  in  S.,  . 

The  proof  depends  on  the  following  lemma,  for  which  we  recall  the  abbreviation 

Qv  s-  (v(1),v(1)+1) *[v(2),v(2)+1) 
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used  in  Section  2  for  the  particular  M*h  square  whose  lower  left  corner  is  the  vertex  v  . 

Ism  4.1.  Let  0  i*  conv{0,  jd^ ,  jd^+dj,  d2>  ,  0  -  U  ®( v,0)  •' 

with  0  the  triangle  conv{jd1#  jd^+dj,  (j+1 )d^ }  .  Let  X  «“ 

{f|Q  (  f  e  st  .  supp  fC  (x(2)  >  l}Vjn  >  .  Then  dim  X  «  (j+1-u>+  . 

Proof.  Since  k  -  3U-2  and  k  -  p  •  H  ,  we  conclude  fro*  (2.5)  that  S(<k, p,  p, p)  is 
spanned  by  the  single  cone  spline  C  :*  C_  with  a  (d^  :  u«  djip,  d^! u)  .  Thus,  from  the 
argument  for  Proposition  2, 

X  -  {f  e  span(c<—  *  f|6“  0  >  • 

The  cone  spline  C  is  homogeneous  of  degree  k  “  3p-2  and  vanishes  to  exact  order 
p  “  2p-2  along  .  Therefore 

C(x  -  (v,0))  -  c  (x  -  (v,0))<M"1'2l*"1)  +  o(x(2)2m_1)  for  xe  0 

for  some  nonzero  c  •  The  condition  *  0  therefore  implies  the  condition 

j  I  ® 

I  a  (x(1 )  -  v)1*"1  »  0  for  x  e  0  (1) 

v-0  v 

in  case  f  -  avC(*-(v,0))  .  Since  the  (univariate)  polynomials  (•»«)**"*,  w0,...,|i-1, 

are  linearly  independent  over  any  open  set,  (1)  constitutes  min {j+1 , p)  independent 
conditions  on  the  coefficient  vector  (ay)  ,  and  therefore 

dim  X  <  (j+1  -  u)+  . 

The  reverse  inequality  follows  from  the  fact  that,  by  Proposition  1 ,  the  box  splines 
M( •  -  ( v,0 ) )  ,  v“0,...,j-j»  ,  are  independent  over  0  and  their  restriction  to  Q  lies 

in  X  .||| 

Corollary.  If  f  «  Sj  has  its  support  in  H,.2  ,  and  its  support  in  [0,H)2  lies 
between  the  rays  (|i-1)d^  +  R+dj  and  ydj  +  R+dj  ,  then  f  vanishes  on  all  of  [0,M)2  . 

Proof.  The  given  domain  lies  in  the  union  of  the  sequence 


of  sets 

Qv  »“  (v,v)  +  conv{0,  (y-Dd^,  (iMld^+dj,  d^} 

8V  (v,v+1)  +  conv{0,  ( M— 1 ) d  ,  (|i-1)d  +d  ,  d  } 
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I 


(■■  illustrated  in  Figure  4.1  for  u  m  2  and  M  ■  4)  to  which  we  can  apply  Ltn  4.1  with 
j  ”  u-1  in  sequence,  in  order  to  conclude,  step  by  step,  that  f  must  also  vanish  on  each 
set  in  the  sequence.  | | | 


Figure  4.1 


4.1, 


Proof  of  Anne  4.1.  Assume  that  supp  g  ^  supp  H  for  some  g  e  . 

with  j  -  u  ,  implies  that  g  »  CM  on  0  i“  conv{(  u-1 ,0 ) ,  (|i,0),  (2|i,U),  (2|i,|it1)}  for 
soom  scalar  c  .  Thus  f  i«  g  -  cm  has  support  in  the  hexagon  supp  M  V.  0  .  This  hexagon 
lies  in  a  domain  of  the  type  described  in  the  corollary  to  henna  4.1,  thus  allowing  the 
conclusion  that  g  ■  cM  .  | | | 


i 

[• 
I 
E 


The  unique  minimality  of  the  support  of  N  and  its  translates  implies  that  they  form 
a  basis  for  the  locally  supported  functions  in  Sj  . 

O 

Proposition  4.1.  For  given  convex  0  ,  {M(  *-v)  :  v  e  Z  ,  supp  M(  *-v)  0}  is  a 

basis  for 

IQ]  :*  (f  C  s  supp  f  C  0)  • 

The  proof  of  this  corollary  is  analogous  to  the  slightly  more  complicated  proof  of 
Proposition  4.2  below  and  is  therefore  omitted. 
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We  now  consider  the  slightly  more  complicated  case  k  -  3p-3  with  p  -  p(k)  ■  2  p-3  . 

Set 

•.  -  «£.,  ■  '!£.«  • 

The  box  splines  in  Sq  do  not  have  minimal  support.  But  an  element  N  of  unique  minimal 


support  in  S0  is  given  by  the  rule 


if  x  e  conv  {0 ,  ,  dj  } 


;  0  otherwise 


Vi  s“  VM1,1,1  for  ve" 


Figure  4.2 


Here,  *  indicates  convolution.  The  function  Ny  seems  to  have  been  considered  first  by 
Frederickson  [Fr]  and  later,  independently,  by  Sabin  [Si] ,  and  thence  in  [SI]  and  [Fa] . 
While  the  support  of  has  seme  symsetry,  it  is  asysnetric  with  respect  to  A  .  Figure 

4.2  shows  the  support  of  Nj  .  In  general,  the  support  of  is  circumambulated  by 

walking  alternatively  p  and  U-1  steps  in  the  directions  d,,  d3,  d2,  -d,,  -dj,  -d2, 
starting  at  the  origin.  Because  of  the  asymmetry,  the  element  given  by  the  rule 

N|/x)  s-  H||(x(2),x<1))  (4 

is  essentially  different  from  .  Together,  they  provide  a  local  basis  since  they  are 
closely  related  to  the  box  splines  in  Sq  j  Convolving  the  obvious  identities 

■*  +  "i  "  *1,1,0  '  ",  +  ■?<-«,>  -  *1,0,1  '  V~V  +  "i  -  *0,1,1 
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with  M 


1,1,1  ' 


obtain  the  identities 


Ny  +  “’m 

"p  +  Vd1> 

VV  +  Nu 


U,  V,  P-1 

*W,  p-1,u 


l**1.  U,  M 


Theorem  4.2,  Let  k 


-  3p-3  and  p  -  p(k)  -  2 p-3  ,  and  let  S0  «■*  *£  ^  •  Then  H ^ 


given  by  (3)  hae  unique  minimal  support  in  sn  . 


The  proof  of  Theorem  4.2  is  based  on  the  following  lemma  and  its  corollaries. 


4.2.  Let  SI  and  X  be  as  in  Lemma  4.1.  but  with  S1  replaced  by  SQ  .  Then 


dim  X  -  (j+1  -  M)+  ♦  ( J+2  -  U>+ 


Proof.  Since  k  -  3y-3  and  k-p  ■  v  ,  we  conclude  from  (2.S)  that  S(<k,p,p,p)  is 
spanned  by  the  two  cone  splines  Cq,  C1  corresponding  to  the  direction  multiplicities 
(u,P-1 ,  M)  and  <W-1,y,U)  ,  respectively.  Thus,  from  the  argument  for  Proposition  2, 

x  -  {f  e  spM(ci<..cv,o),|0)i:0fVi,  «  f|e-  0}  , 

with  0  :»  conv{ jd^  ,  jd^+dj,  ( j+1 ) d^ )  .  Because  is  hasK>geneous  of  degree  3j*-3  and 

vanishes  to  exact  order  p+i  “  2y-3  +  i  along  Ml  ,  we  have 

C0<x-(v,0))  -  c0(x-(v,0))<l>■1'2,,"2,  +  d0(jv-(v,0))(,l"2'2,,"1>  +  o(x<2)2M"1  ) 

C,  (x-(  v,0>)  -  c1(x-<V,0))<‘,■2,2,,"1,  +  o(x(2)2M_1) 

for  x  e  0  and  some  cQ,  c,  j*  0  .  The  condition 


f  i-  E  a  -  <v,0))  -  0  on  0 


therefore  implies  that 


E  «  0(x(  1 )  -  v)» 


for  x  e  0 


These  ere  min (j+1 ,  p}  +  ein{ j+1 , y-1 }  lineerly  independent  conditions  on  the  2(  j+1 )  - 
vector  (a^)  of  coefficients  end  therefore 

di*  X  <  (j+1  -  M)+  +  ( j+2  -  u)+  . 

The  reverse  inequelity  follows  since 

*V*  ”  •  WO . j-**,  end  •  -  (0,v))  ,  \M), . . . ,  j+1-u 

ere  independent  over  Q  end  their  restriction  to  Q  lies  in  X  .  Their  lineer 
independence  follows  from  the  feet  thet,  by  (5),  end  agree  near  the  origin  with 

CQ  and  C1  ,  respectively.  | | | 

Corollary  1.  Let  0  end  0  be  as  in  the  lata,  but  with  j  ”  u-1  ,  end  let  O'  end 
0'  be  their  image  under  the  'flip'  x  | — >  (x(2),x(1))  .if  f  e  SQ  has  support  in 
{x(1),  x(2)  >  1>UQU0a  ,  then  f  vanishes  already  on  Qv/Q*  . 

Proof.  Much  as  in  the  proof  of  Leona  4.2.  we  conclude  that,  on  (IV/  Q'  , 

f  “  a00C0+a01C1  +  E  (aviCi(-(v'0>>  +  aUci(-(0'v,)) 

0<  v<(i 

i«0,1 


subject  to  the  conditions  that 


E  a  (x(1)  -  v)**”1  ■  0  for  x  e  0 


1  (\odo  +  aviciKx<1)  ~  V^2  ‘  0  *  e  8  ( 

v<u  ™ 

for  sone  c^  /  0  ,  while,  with  aj^  «-  agA  ,  1-0,1,  also 

E  *’  (x(2)  -  *  0  for  x  e  O'  (6 

v<y 

1  (*Cld0  *  a^0clKx<2>  "  v)  -  0  for  x  e  0a  .  (6 

V<JI 

Dote  the  reversal  in  the  role  of  the  second  subscript,  due  to  the  fact  that 

C^x)  -  C1_i(x(2),x(1 ))  . 

He  conclude  froa  (6a)  that  a^  -  0  ,  all  v  ,  and  frosi  (6'a)  that  a^  -  0  ,  all  v  .  In 
particular,  a£,  -  aQ1  -  0  .  Therefore  < 6b)  ioplies  that  also  ay1  -  0  ,  all  v  ,  and, 
likewise,  (6'b)  ioplies  that  a'  -  0  ,  all  v  .  ||| 

VO 


29 


sJ 


OwolUry  2.  If  f  (  Sg  has  its  support  In  ^2  and  its  support  In  [0,M]2  lias 
between  tha  rays  (u-Dd^  +  R+dj  and  (i*-1)d2  +  R+dj  ,  then  f  -  0  on  [0,Ml2  . 


Proof.  Tha  gi van  domain  ia  eontainad  in  tha  union  of  tha  aaquanca 


of  sata 


V  °1*  ***'  °H-1 


oy  «-  <v,w)  +  ava*  , 


with  alia*  aa  in  Corollary  1,  to  which  we  can  apply  Corollary  1  in  aaquanca,  in  order  to 
conclude  that  f  auat  vanish  on  each  .  | | | 


Proof  of 


4.2.  Assume  that  aupp  g  £  aupp  n  for 


9  €  Sq  •  By  Imt 


4.2,  g-  cN^  on  8  i»  conv{ (m-1  ,0 ) ,  (p,0),  (2p-1,p-1),  (2p-1 , M> >  .  This  implies  that  f  j* 
g  -  CK  e  Sq  has  support  in  tha  doauin  described  in  Corollary  2  to  Lenma  4.2,  hence  must 
be  aero.  | | ] 


we  now  show  that  tha  elements  with  unique  minimal  support  form  a  basis  for  all  locally 
supported  elements  of  SQ  . 


Propoaitiom  4.2.  For  given  convex  8  ,  a  basis  for  SqIOI  i*>  {f  e  S0  :  aupp  f(0)  is 
provided  by  tha  collection  of  all  N  («-v),  W  <  »-v)  ,  v  e  C2,  in  s  [8]  . 


Proof.  Assume  without  loss  that  A  £  [0,M]  and  let  f  e  SQ[8]  be  given.  Since  , 
PJ,  agree  near  tha  origin  with  tha  respective  cone  splines  CQ,  C,  ,  tha  argument  for 
Proposition  2  leads  to  the  conclusion  that  thars  is  a  unique  linear  combination 


veiO,M-2t*1]  9  9 


which  agrees  with  f  on  [0,M-2i*+2]2  .  This  implies  that  f-g  has  support  only  on 
[0,M]  2\  (0,M-2li+-2]2  .  Application  of  Corollary  2  to  I  earns  4.2  therefore  proves  that  f-g 
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has  no  support  in  {x(1)  <  m-2ji+2>  and,  with  thia,  a  second  application  of  that  corollary 
shows  that  f-g  has  no  support  in  {M-21H-2  <  x(1)  <  m}  either.  Therefore  f  -  g  • 

Let  Q'  be  the  convex  hull  of  the  union  of  the  supports  of  all  the  M^( *-v)  and 
Nj|( -v)  which  appear  in  (7)  with  nonzero  coefficients.  Then  S'  is  a  polygon.  We  claim 
that  0*  C  8  .  It  is  obvious  from  the  construction  of  g  that  any  lower  left  corner  of 
8'  must  lie  in  supp  f  ,  hence  in  8  .  But  since  g  is  uniquely  determined,  this 
implies,  using  the  mesh  symmetries,  that  all  six  kinds  of  corners  of  8*  lie  in  8  ,  hence 
so  does  8'  .  Ill 


Finally,  we  consider  the  irregular  and  rich  case  k  “  3u-1  ,  for  which  p(k)  “  2 p-2  . 


Set 

S2  k,A  3u-1,A  * 

There  are  three  independent  cone  splines  of  degree  k  at  a  lower  left  corner,  but,  because 
p(k)  is  so  low,  there  is  also  an  additional  cone  spline  of  degree  k-1  in  S2  •  This  means 
that  a  search  for  a  basis  for  S2 (8)  would  have  to  come  19  with  four  unique  minimal 
support  elements  per  vertex.  This,  as  it  turns  out,  is  not  possible  if  we  stick  with  the 
definition  of  "minimal  support"  given  earlier. 


Already  the  case  k  -  2  of  continuous  parabolic  splines  provides  the  necessary 
illustrations t  A  suitable  basis  for  S2 [Q]  is  provided  by  the  translates  of  four  functions 
whose  supports  are  drawn  in  Figure  4.3.  These  functions  are  obtainable  as  the  Lagrange 
functions  of  standard  local  parabolic  interpolation  (at  the  vertices  and  the  edge  midpoints 
of  each  triangle).  The  first  three  have  unique  minimal  support.  But  the  fourth  function's 
support  is  made  up  of  six  triangles  and  could  accomodate  each  of  the  other  three's  much 
smaller  support,  hence  it  fails  to  be  minimal. 

In  this  case  and  others  mentioned  later,  it  is  possible  to  recapture  unique  minimality 
of  support  by  referring  to  the  support  of  the  associated  Bternstein  or  -ezier)-net  of  the 
pp  functions  instead.  In  any  case,  the  support  of  the  various  box  splines  of  degree  2  in 
S2  is  far  from  minimal  since  it  contains  ten  triangles. 

The  next  case,  k  -  5  ,  hence  p(k)  -  2,  provides  the  additional  unhappy  surprise  that, 
in  this  case,  S2  contains  an  element  supported  on  just  one  hexagon  (i.e.,  on  six 
triangles).  This  element  occurs  already  in  [SI].  Because  its  support  coincides  with  that 
of  M1  j  ^  ,  this  element  cannot  be  obtained  from  a  parabolic  one  by  convolution.  The  same 
is  true  of  the  next  two  "minimal  support"  elements  whose  supports  coincide  with  that  of 
and  respectively  (see  Figure  4.2).  The  fourth  "minimal  support"  element  is 

derived  from,  and  has  the  same  support  as,  Mj  2  2  . 

The  next  case,  k  «  8,  is  easy  since  its  four  "minimal  support"  elements  can  be 
obtained  from  those  for  k  «  5  by  convolution  with  . 

This  pattern  repeats:  For  odd  u  ,  the  four  "minimal  support"  elements  can  be  obtained 
from  the  preceding  case  by  convolution  with  .  For  even  u  ,  enough  local  degrees  of 

freedom  have  been  accumulated  to  make  possible  elements  in  S2  of  yet  smaller  support  than 
is  had  by  the  elements  obtained  from  the  preceding  case  by  convolution  with  M.  ,  .  . 

'MM 
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